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Abstract 

In this paper, we consider the long time asymptotics of multi-time correlation functions 
for quantum dynamical systems that are sufficiently random to relax to a "reference state" . 
In particular, the evolution of such systems must have a continuous spectrum. Special 
attention is paid to general dynamical clustering conditions and their consequences for the 
structure of fluctuations of temporal averages. 
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1 Introduction 



One of the main problems in quantum chaos is to understand the relaxation phenomena 
induced by the dynamics in systems with few degrees of freedom. Typically, the system 
relaxes on an appropriate time scale characteristic of the dynamics. Indeed, the spectrum of 
the evolution of chaotic quantum systems is usually discrete and observation of the system 
for very long times will reveal this discrete nature: the time correlation functions are quasi- 
periodic. The separation between the points of the spectrum depends on a quantization 
parameter such as Ti in the Chirikov kicked rotator [1] or on the dimension of an irreducible 
representation of SU(2) such as in the kicked top [2] or on that of an underlying Hilbert 
space as in the finite-dimensional Cat maps [3, 4]. When the quantization parameter tends 
to an appropriate limit, one obtains a classical dynamical system. We are not concerned 
with the statistical properties of eigenvalues and eigenvectors of such models, but rather 
with dynamical properties of expectations values. It is a typical feature that the classical 
limit and the limit for large times of such expectation values cannot be exchanged and so 
one can look for a joint limit obtained by a suitable scaling of the time with respect to the 
quantization parameter [5]. The aim is to extract by such a procedure true relaxation in 
rescaled time. 

We don't address this scaling problem in this paper, but rather concentrate on the large 
time behaviour of model systems with already fully displayed relaxation as a consequence 
of basic dynamical randomness. Different types of behaviour are possible and we introduce 
a way of describing them by considering the asymptotic analysis of multi-time correlation 
functions [6]. In particular it appears that fluctuations around temporal averages are 
a useful tool to distinguish between various degrees of randomness. For instance, the 
distribution of fluctuations may be given by the usual Gaussian law, but also by the 
semicircle law or by more exotic laws. Such distributions have already been obtained 
considering fluctuations of observables obeying stochastic commutation relations [7]. In 
this paper, we show that deterministic dynamics, i.e. without any stochastic input, can lead 
to quite a variety of distributions for fluctuations. The appearance of unusual statistics, 
such as the free statistics associated with the semicircle law, is connected with chaotic 
features of the dynamics. By chaotic or random quantum features, we shall mean different 
degrees of clustering in time, namely different strengths with which events largely separated 
in time tend to become independent. Classically, one has the notion of mixing, whereas 
in quantum mechanics two notions are commonly considered, that of weak and strong 
clustering, the latter one implying the former. It turns out that the appearance of exotic 
statistics at the asymptotic level, rater than the common, Gaussian one, associated with 
the notion of classical independence and arising from the stronger clustering properties, is 
related to a finer distinction of possibilities between weak and strong clustering. We anyway 
expect to observe the emergence of such statistics, not only in dynamical systems with fully 
displayed relaxation, but also for fluctuations in appropriately scaled finite systems. 

In Section 2, we review some notions of randomness for quantum dynamical systems. 
Section 3 deals with the construction of the asymptotics of multi-time correlation functions 
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and provides hereby a useful setting for describing the law of large numbers for a sufficiently 
ergodic dynamical system. Then, we consider in Section 4 how a central limit theorem 
for fluctuations can be obtained. In Sections 5, we present the structures that arise when 
considering some simple examples of dynamical systems. 

In this paper, we shall model quantum dynamical systems by triples (21, 0, 4>) where 

• 21 is a unital C*-algebra of observables 

• © = {©t t £ R or Z} is a dynamical group of *-automorphisms of 21 either in 
continuous or in discrete time and X(t) will denote the observable X e 21 evolved up 
to time t: X(t) = Q t {X), and 

• the reference state is assumed to be invariant under 6 i.e. <j) o Q t = (p. 
General multi-time correlation functions are functions of the form 

t^^iW(^ (1) )i( 2 )(f, (2) )...i(")(t, (fl) )) , (1) 

where the X^\t u ^) are operators at times t v ^ in 21, t = {£i,t2, • • •} G Z N ° and v maps 
{1,2, . . .,71} into N . 

Usually, in quantum statistical mechanics, one considers time-ordered correlation func- 
tions. Since one expects observables largely separated in time to commute, such an ordering 
would be no restriction. We are interested, however, in a situation where a complex dy- 
namics generates wildly fluctuating algebraic relations for observables largely separated 
in time. In such a case, commutation relations cannot be used to simplify correlation 
functions by grouping together observables at equal large times and, in order to obtain 
information about algebraic relations between observables largely separated in time, we 
have to consider general expressions as in (1). In particular, it is necessary to include the 
possibility of repeating a same time label within a correlation function. Then, the natural 
algebraic structure to consider is that of a countable free product 21^ of copies of 21 [8], 
while the asymptotics of multi-time correlation functions will allow us to equip 21^ with 
an asymptotic state 0oo. The probabilistic structure corresponding to such a state should 
reflect the essential features of the underlying dynamics. 

We briefly remind the construction of 21^ as the free product *ieN 2li and we shall 
refer to 21^ as the asymptotic free algebra associated with the dynamical system (21, 0, 4>). 
Each of the algebras 2tj is a copy of the basic algebra 21 of observables and 21^ is the univer- 
sal C*-algebra generated by an identity element 1 and by "words" w = X^X^) • • • X^^ 

that consist of concatenations of "letters" X^ e 21. The subscript v(j) in X^^ refers 

to which copy of 21 the letter X^ belongs to and any two consecutive subscripts are 
unequal. Concatenation, together with simplification rules, defines the product of words. 
More specifically, the rules for handling words are: for X, Y e 21, A G C, j G No and w, w' 
two generic words 

wljw' =^ ww' (2. a) 
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w (Xj + XYj)w' wXjw' + XwYjw' (2.b) 
wXjYjw' => w(XY) jW ' . (2.c) 



Notice that the product XY in (2.c) is not concatenation, but rather the usual operator 
product in the algebra 21. Moreover, the adjoint w* of a word w = X^^X^^ ■ ■ ■ X^^ 

equals (l W t w (X M t M -(X (1) t (1) . 



2 Random behaviours in quantum systems 

Before considering the problem of endowing St^ with an asymptotic state, we present 
a hierarchy of ergodic properties typical for infinite quantum systems [9, 10]. We shall 
formulate them for the case of discrete time dynamical systems (21, ©, 4>). 

Actually, our construction will be essentially based on the use of time averages of 
correlation functions of the form 

<P(XY(t)Z) aV := lim i ]T <j>{X Y{s) Z) , (3) 

where X, Y and Z are arbitrary observables in 21. 

If the dynamics of a system is sufficiently regular, observations turn out to be quite 
correlated, even when largely separated in time. We shall then associate various degrees 
of randomness with the strength of decorrelation properties of the dynamics, if any. The 
lowest degree of randomness is clustering in the mean 



(f>(XY(t)Z) =<f>(XZ)(f>(Y), X,Y,Ze%, (4) 
next in the list comes weak clustering 

\mi(j>(XY(t)Z) = (j>(XZ)(j>(Y), X,Y,Ze%, (5) 

and we close the list with strong clustering 

lim <f)(XY(t) Z S(t)T) =<f)(X ZT)<f)(Y S), X,Y,Z,S,Te% (6.a) 

and hyper-clustering 

in[|ti ^,_ h* ,,, <w* < *><^>) ■ -n*flL-. w * , " ) ) ■ (6 - b) 

3 

In formula (6.b), the limit is taken in such a way that all the times t v (j) and the differences 
between those with different indices go to infinity. Notice that, as in (1), repeated times 
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are allowed in (6.b) and the arrow over the product at the right-hand side indicates that 
the order among operators at equal times has to be preserved since they do not commute 
in general. Among the previous properties, one can check that the following relations hold 

(6.b) <S> (6.a) (5) =>■ (4) . (7) 

While the implications from left to right are evident, in order to deduce the equivalence 
in the first place, we use the notion of asymptotic Abelianess in time. Indeed, the different 
kinds of clustering (4), (5) and (6. a), imply the following degrees of asymptotic Abelianess: 
asymptotic Abelianess in the mean 



(4) =^ (j>(S [X, Y(t)] Z) =0 VS, X, Y, Z G 21 , (8.a) 
weak asymptotic Abelianess 

(5) =>. lim <f)(S [X, Y(t)] Z) = VS, X,Y, Z G 21 (8.b) 

and strong asymptotic Abelianess 

(6.a) =^ lim (j)(S* [X,Y(t)]*[X,Y(t)]S) =0 VS, X, Y G 21 . (8.c) 

t— >oo 

Weak clustering plus strong asymptotic Abelianess imply strong clustering. Indeed 
lim cj)(X Y(t) Z S(t) T) = lim <f>(X Y(t) S(t) ZT)+ lim <f>(X Y(t) [Z, S(t)] T) 

t — >oo t — >oo t — >oo 

= <j)(X ZT)cf)(Y S) VS, T, X,Y, Z G 21 . 

The second limit tends to zero because of strong asymptotic Abelianess and the Cauchy- 
Schwartz inequality 

\(j)(XY(t)[Z, S(t)} T)\ 2 < (j)(XY(t) Y(t)*X*) (j)(T*[Z, S(t)]*[Z, S(t)} T) 

< \\X\\ 2 \\Y\\ 2 (P(T*[Z,S(t)]*[Z,S(t)]T) , 

while the first one gives the result because of weak-clustering. By a similar argument, (6. a) 
implies (6.b) and therefore the equivalence in (7) holds. We just sketch here the main idea 
of the argument by considering the case of three equal times. In a first step, for fixed 
X, Y, Q, S, T, U and Z in 21, write 

<f>{X Y(t) S Z(t) T U(t) Q) = (p(X S Y(t) Z(t) T U(t) Q)+ 

cj>{X[Y{t),S]Z{t)TU{t)Q) . (9) 

By the Cauchy-Schwarz inequality and (p(abb*a*) < \\b\\ 2 <f>(aa*) we obtain 

\4>{X [Y(t), S] Z(t)TU(t) Q)\ 2 < \\Z\\ 2 \\T\\ 2 \\U\\ 2 \\Q\\ 2 cf>(X [Y(t), S][Y(t), S]*X*) . 

Assuming strong asymptotic Abelianess (8.c), the second term in (9) vanishes when t 
tends to infinity. Repeating the same procedure on the first term, we collect together 
(YZU)(t) and the conclusion follows from weak-clustering. Indeed, the condition that 
inf \ti — tj\ — > oo means that all different times are so largely separated that the regrouped 
correlation functions cluster asymptotically with respect to any of them. 
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3 The asymptotic multi-time correlations 

Consider, as in (1), a multi-time correlation function 

tH^I«(t Kl) )l( 2 )(t K2) )...lW(^))) . (10) 

We may restrict the map v : {1,2, ... ,n} — > No in such a way that if a u(j) appears 
somewhere in the operator product in (10), then all smaller natural numbers must have 
already appeared at least once to the left of v(j). More formally 

if 1 < i < then there exists 1 < i < j such that v(i) = £ . (11) 

As a consequence of this prescription, z/({l, 2, . . . , n}) = {1,2, ... ,k} with k < n due 
to possible repetitions of a label. Notice that this is a convenient way of rewriting the 
multi-time correlation functions that were introduced in (1): e.g. t i— > c[)(X (t^)Y fa)) is 
not excluded by the rule of above for v since it equals t i— > (p(l(ti)l(t2))X(t3)Yfa)) . 
The reason for making this specific choice in writing correlation functions is to have all 
labels 1,2, ...,n appear, in increasing order and possibly with repetitions, in a correla- 
tion function if the label n appears. This will become useful later on when we consider 
time averages. The multi-time correlation functions (10) form a self-adjoint linear space 
containing the constant function. 

We shall describe the asymptotics of multi-time correlation functions in terms of a 
state (^oo on the asymptotic free algebra 2loo. tp^ is obtained by averaging multi-time 
correlation functions. We now generalize the procedure of averaging considered in (3) for a 
single-time correlation function and briefly recall the notion of invariant mean [11]. A mean 
Avg on a set X is a normalized and positive linear functional on a self-adjoint linear space 
X of bounded complex-valued functions on X containing the constant function and closed 
for the || Hoo-norm. If is a group of transformations of X, leaving X globally invariant, 
then Avg is called invariant if Avg(/) = Avg(/ o 7). 

As we consider discrete time dynamical systems (21, O, 4>), we take Z N ° for the space 
X . The group (25 will consist of multi-time translations on X given by 

7s :A^A:t^t-s, s = {s x , s 2 , . . .} G Z N ° . 

The function space X on which we shall consider means is the closure of the linear space 
of multi-time correlation functions. X is invariant under the group action / 1— > / o ^y s since 
is an automorphism of 21. 

In general, we want to explore the possibility of defining a positive, normalized linear 
functional (a state) (j)^ on the linear span of the words in the asymptotic free algebra, via 
some mean "Avg " defined on the space of multi-time correlation functions. 
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We proceed as follows: first we rewrite a general word w = X,^X^ 2 X. • ■ ■ X,, 
without restrictions on the map \x except that consecutive subindices \x{j) and + 1) 
are different. By inserting an appropriate number of identities we can write in a unique 

way w = 

r ia) y i 2 2) • • - y !aIv where the Y 

are either identities or are X's and where v 
satisfies (11). E.g. w = X2 is rewritten as w = I1X2 and w = X2Y1Z2W 4 becomes 
w = I1X2Y1Z2I3W4. Given the space of multi-time correlation functions of a dynamical 
system (21, 6, (p), we shall consider, when existing, the linear functionals 

0oo(4ji)4S)'--4S)) :=Avg(t^^(lW(t Kl) )l( 2 )(f K2) )...lW(^)))) , (12) 

where v maps {1,2,..., n} into {1,2,..., s}. Furthermore, we restrict ourselves to those 
means Avg on X that satisfy the following strong compatibility condition: 

if (£1, t2, . . . , t n ) I— > f(ti, t2, . . . , t n ) is a uniform limit of multi-time correlation 
functions, and if / depends only on a subset {£i(i), £i(2)> • • • > ^i(fe)} °f variables 
where i : {1,2, ... ,k} — > {1,2, ...,n} is an order preserving injection, that is 
f(t 1 ,t 2 ,...,t n ) =g(t i ( 1) ,t i{ 2),...,t i{k) ), then Avg(/) = Avg(^). (13) 

Proposition 1. if Avg existe as an invariant mean on the space of multi-time correlation 
functions and satisfies the strong compatibility condition (13), then the functional (j)^ 
defined in (12) extends to a state on SIqo. Furthermore, 



o 



S = (p^ and 4>oo do = (poo , 



where S satisfies S (X,) := (0 Sj . (X))_. for j G N and X e 21, whereas 6 is any order 
preserving infective transformation o/N andag is the *-homomorphism o/SIqo determined 
by a e (Xj) := X e ^y 

Proof: Normalization and linearity are a consequence of (2. a) and (2.b). In order to prove 
positivity we have to check that for any finite linear combination W = J2 W X w w of words 
one has (poo(W*W) > 0. From definition (12) it turns out that (poo(W*W) is the multiple 
average of the expectation of a positive operator in the state (p, hence positive. Multi-time 
invar iance follows from the invar iance of the mean. Again, instead of producing a formal 
proof of the invariance of (poo under ag, we present a simple example which clarifies the 
essential mechanism. We show that (^(A^Yi-^) = (pooiX^YiZ^), in which case 9(1) = 1 
and 9(2) = 3. Relabelling dummy variables and using condition (13), we obtain 

tooiXsYtZs) = Avg((t 1 ,t 2 ,t 3 ,...)^ ( p(l(t 1 )l(t2)X(t 3 )Y(t 1 )Z(t 3 : 

= Avg((*i, £2, £3, • • •) - 0(l(£i)X(£ 2 )n*i)^2))) 
= ( p oo (X 2 Y 1 Z2) . 
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In the following, we shall fix a definite averaging procedure and construct asymptotic 
states averaging over the different times in successive order 

Avg(t ^ /(*l,*2,---,*n)) : = 

Avg(t n i-> •••Avg^ >-> Avg^i h-> /(*i,t 2 , - --,*«))) •••) • (14) 

More general procedures are possible by considering coupled limits: e.g. one could try to 
construct a mean by averaging functions (t\, £2, • • • , t n ) 1— > f(ti,t2, ■ ■ ■ , t n ) of n arguments 
over n-dimensional cubes 

^ L L 

Avg^(ti,t 2 ,...,t n ) ^ f(ti,t 2 ,...,t n )^) := 2 j™ ) ^ "' /(*!>•••>*») • 

ti=i t„=i 

Such a mean is also strongly compatible. Notice, however, that, choosing different single- 
time averages Avg in (14), we can in general violate the strong compatibility condition (13). 
Here follows an example with two times and two different single-time averages 

Avg 2 (t 2 >-> Av gl (ti i-> /(£ 2 ))) = Avg 2 (t i-> /(£)) 

^Avgi (*->/(*)) 

= Avg 2 (t 2 i-> Avgi^i i-> /(ti))) • 

Quite to the other extreme, there are examples of dynamical systems and of averages 
of their multi-time correlation functions that are not only strongly compatible but even 
permutation invariant in the sense that 

Avg(a w (/)) = Avg(/) (15) 

for any local permutation n of the natural numbers, that is for any bijection 7r from No 
into No that leaves all elements invariant except for a finite number of them. a n acts 
on a correlation function t 1— > /(£i,t 2 ,...) as ct n (f)(ti, t 2 , . . .) = /(t ir (i),t 7r (2)> • • •)• As a 
consequence of (15), the asymptotic states (j)^ defined by permutation invariant multi-time 
averages will also be permutation invariant on the asymptotic algebra 2too, namely 

<j> 00 (xW x w...x ( "\) =<j )oo (x il) ( n X (2) (0 ,---X {n) ( ,) VX^eSt. (16) 

V i/(2) i/(n) 7 V nou(l) 7roi/(2) 7ro;/(n) / V / 

The asymptotic states 0oo are always permutation invariant when the basic dynam- 
ical system (21, 6, (j>) is strongly clustering. In order to prove this, we introduce a useful 
technical result. 
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Lemma 1. For d,k G N define 

A h, / \ f zf Iti — t»| < d for some 1 < % ^ j < k 

A d (t 1 ,...,t k ) = { i 2 e jl - - rj - ■ 

Then, if the multiple average (14) of a uniformly bounded function f : Z k — > C exists, we 
have 

Avg(t k ^■■■Avg(t 2 i-> Avg^i i-> /(*i,t 2 , - - -,tfc))) •••) = 

Avg^ fc i-> •••Avg^ h-> Avg^i ^ A3(ti,...,tfc)/(ti,t2,-..,tfc))) •••) • 



Proof: For fixed d, choose t±, . . . , G Z in such a way that the function 

A^ _1 (ti, . . . .,t fc ) = 1. Then, A5(*i, . . . , . . . , t k ) = for only a finite num- 

ber of values of U G Z which implies that 



Avg^j i-> A5(ti,...,ti, ...,tfc)^(ti, ...,tfc) 

= Aj -1 (*i, . . . , *<_i, tt+i, • • • , tfc) Avg (ti i-> g(ti, ...,U,..., tkfj 

for a uniformly bounded g : Z fc — > C. By successively applying this observation to 

Avg(t n i-> ■■■Avg(t 2 h-> Avg(ti h-> A3(ti,..., t n )/(ti, t 2 ,...,*»))) •••) , 

the lemma follows. Whereas in the statement of the lemma a definite order of the time 
averages has been specified, its proof is independent of it. 



Proposition 2. Let (21,0,0) be strongly clustering as in (6. a), equivalently (6.b). Then, 
if the multiple average in (14) exists, it defines a permutation invariant asymptotic state 
4>oo on 2loo . 

Proof: Let us consider a multi-time correlation function 

cj>[x^\t v{1) )X^\t v{2) ) ■ ■ ■ X^\t v{n) )) , 

where, because of strong compatibility, we can assume that v : {1, 2, . . . , n} — > {1, 2, . . . , k} 
with k < n. Thus, we collect all i G {1, 2, . . . , n} such that v(i) = p into subsets I p . Because 
of strong-clustering, and thus of hyper-clustering (6.b), for any e > 0, we can choose d > 
such that \ti — tj\ > d for all ti ^ tj, with i,j G {1, 2, ... , k}, implies 



K > 

^(i (i) (Mi))i (2) (^)) • ^ (b) m) - n <Kn 67 x(i) 



P =i 



< e . 
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Thus, when evaluating p^i^X^^X^y ■ ■ X^ n ^ via the prescription (14), we can use 
Lemma 1 and a function A^(ti, £2, • • • , tk) to estimate 



'°°V ^-1/(1)^-1/(2) ' ' i/(n) 



n*0L.* w 



p=l 



< e . 



Therefore, in the case of strong clustering, in whichever order the single time averages 
are performed, the multiple average (14) agrees with the time limit (6.b) of multi-time 
correlation functions. The ensuing asymptotic state p^ is thus permutation invariant. 
Moreover, from time invariance of p, we have that p 00 (Xj) = p(X) for any j e N and 
X E 21. Therefore, with i e I p implying v(i) = p, 

^(xM ) ...x<») ) )=n^(n jel x(<)). (17) 

P =i 



When the basic dynamics is strongly clustering, the structure of the expectations on 
2loo calculated with respect to the asymptotic state p^ corresponds to the usual notion 
of commutative independence of random variables. Indeed, p^ vanishes on the two-sided 
ideal of 2loo generated by commutators of letters sitting in different copies of 21 in 2loo, i.e. 
by {LYj,Yfc] \ j ^ k, X,Y e 21}. We may therefore think of p^ as an infinite product 
state on the minimal tensor product of a countable number of copies of 21. 

A very different notion of independence, called freeness or free independence, has been 
recently introduced in the realm of non-commutative probability [8]. This notion corre- 
sponds to the following structure for the correlation functions of a state ^ona free product 
of C*-algebras 93 j 

*(4»*g>... *£>)=o 

whenever ip(yX^ j = and jk 7^ jk+i f° r au k. Accordingly, one refers to ip as to the free 
product of the states ipj, where ipj is the restriction of tp to 93j. 

Freeness is totally incompatible with statistical independence as in (17). In fact, if 
poo satisfied both freeness and commutative independence, then we would have for any 
centred observable X, i.e. <p(X) = 0, that 

= p oo (x*x*x x 1 ) = P(X*X) 2 . 
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4 Temporal fluctuations 



The two cases of commutative and free independence, presented in the previous section, 
are somehow extreme. Many other possibilities for the structure of 0oo may arise. In 
general, we cannot hope to obtain a comprehensive description of these structures. A 
more manageable framework is provided by "fluctuations". They are natural objects to 
consider as we may think of the state (p^ as determining the Law of Large Numbers for 
the dynamical system (21, ©,0). Fluctuations are thus on the level of the Central Limit 
Theorem. 

Definition. Let N e N and X e 21. A local fluctuation F^{X) is the following element 

0/2loo 

1 N 

F N (X) ^—^(Xt- 0(X)l) . 

V i=l 

It is our aim to discuss the limiting behaviour of Fn(X) when iV tends to infinity, 
namely to study limits of correlations such as 

lim ( p oo (F N (xW)F N (X^)---F N (X^)) , (18) 
and to establish an algebraic central limit theorem. 

In a stronger sense, one may try to reconstruct possible algebraic relations of global 
fluctuations F(X), if any, by means of the functional $ on global fluctuations defined by 

$(F(lW)F(^ (2) )--f(^ (r) )) := lim 0oo (f n (X^)F n (X^) ■ ■ ■ F N (X^f) . (19) 

Notice that, in the case of commutative, respectively free independence, the linear func- 
tional $ on the r.h.s. of (19) is in fact a well-defined state on the algebra of fluctua- 
tions [12, 7] such that they become Gaussianly, respectively semicircularly, distributed, 
non-commutative, random variables. 

We shall restrict our considerations to averages of multi-time correlation functions of 
the form (14) and study in some generality the limit (18) by adapting an argument in [13]. 
We show that a clustering condition stronger than weak clustering (5), but weaker than 
strong clustering (6. a), is sufficient to ensure that only moments of even order contribute 
to the limit joint distribution of fluctuations. 

Proposition 3. Let us assume that the quantum dynamical system (21, G, <fi) satisfies the 
following cluster condition 

lim ^(Z^)^)) • • ' • ^ (n) (^(n))) = , (20) 

int \ti—tj |— >oo V / 
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whenever Y is centred and v maps {1,2,..., n} into N , where, as in (6.b), lim inf \ t ._ t 
means that all times and differences of different times go to infinity. If the multiple av- 
erage (14) exists, the state (j)^ it defines on the asymptotic algebra 2loo is such that, with 
. . . ,X^ in 21 centred observables, 

r = 2n + 1 

(X {1) ...X {2n) ) r-2n 

v 

Y^v^ means that we have to sum over all partitions v of {1, 2, ... , 2n} into pairs ((«i, fix), 
(«2,/32) ; • • • , (ot n , Pn)) 2.e. we choose sites aj < (3j such that v{a.j) = v((3j) = j with j 
running from 1 to n. 



lim <j>jF N {X^)F N {xM)...F N {X^)\ = 



Proof: We have to compute the limit for large N of 

<I>oo(f n (xW)f n (xM)...f n (x^)) = 1 ^ Yl <m< ) ---4i ) ) > ( 21 ) 

ke{i,...,iv}- 

where k = {ki, . . . , k r }. We first concentrate on the contributions to (21) where at least 
one of the indices k±,k2, ■ ■ • ,k r appears only once and show that all of them vanish. Let 
k p be such an index. Due to strong compatibility, we may always assume that the map 
i 1— > ki transform {1,2,..., r} into {1,2,..., s}, with s < r due to possible repetitions of 
an index ki. However, by hypothesis, the index k p appears just once, meaning that the 
time tk p does not explicitly appear in the words 

X^{t kl )X^{t k2 ) ■ ■■X^\t kp _ 1 ) and X^\t kp+1 )X^\t kp+2 ) ■ ■ ■ X^\t kr ) . 

Let X( pS> be a centred observable and use that (p^ is defined by (14) and that the dynamical 
system (21, 0, </>) satisfies condition (20). Then, Lemma 1 guarantees that, given any e > 0, 
we can find a corridor function A^(ti, ■ ■ ■ , t s ) with d so large that 

<p°° Up ■ • •4 p - 1) 4 p) 4 p+1) • • -4 r) ) = 

roo ^ kl fcp_i k p k p+1 k r J 

Avg(t s ^ ■ ■ ■ A\g(t 2 I ^ Avg^i 

<p(x^(t kl ) ■ ..X^\t kp ^)X^\t kp )X^\t kp+1 ) ■ ■ ■X^\t kr )))) ■ ■ ■) 



< e . 



Therefore, the only non-zero contributions come from terms where none of the subindices 
kj appears alone. 

Next, we prove that, in the limit of large N, all those contributions vanish which come 
from terms where all indices appear twice and at least one thrice. As (j)^ is a state, we 
have the a priori estimate 



i^«...x<:>)i<niixM| 
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Let us partition {1,2,..., r} into s groups of equal indices and let r\, r%, . . . , r s denote the 
number of elements in each such group. By assumption, Tj > 2, at least one of the rj > 3 
and T\ + ri + ■ ■ ■ + r s = r > 2s. The contribution of all these terms can be estimated 
as follows: any partition k : {1,2, ...,r} i— > {ki,k%, . . . ,k 8 }, where kj G {1, 2, . . . , N}, 
can be written as k = 9 o v, with v : {1, 2, . . . , r} i— > {1,2, ... ,s} a partition and 9 : 
{1,2, ... ,s} I— > {1,2,..., N} an order preserving injection. There are (^) of such order 
preserving injections. Thus, the contribution of the type considered is bounded from above 
by 

A s is the number of partitions v : {1,2, ... ,r} i— > {1,2, ... ,s} such that each 1 < j < s 
appears at least twice. Such an upper bound tends to zero when N tends to infinity. 

Therefore, we remain with contributions of the form (f)^ (x\^ • • • X^^j where the kj run 
from 1 to AT and each of them appears exactly twice, whence r has to be even, say r = 
2n. These contributions are given by 0oo (^(i) • ' • ^(2n)) ' wnere v '• {!>••• 5 2n} i— > 
{1, . . . , n} is a (not necessarily ordered) pair partition, that is, for any j e {1, 2, . . . , n}, 
there exist two, and only two, indices (aj, (3j) G {1,2,..., 2n} such that v{a.j) = = 
j. In fact, as before: any pair partition k : {1,2, ...,2n} i— > {k\, ki, ■ ■ ■ , k n }, where 
kj G {1,2,..., iV}, can be written as k = 9 o v with v : {1,2,..., 2n} i— > {1, 2, . . . , n} a 
pair partition and 9 : {1,2, ... ,n} i— > {1,2,..., iV} an order preserving injection. There 
being (^) of such injections, the sum of contributions from generic pair partitions k : 
{1,2,..., 2n} i— > {1, 2, ... , iV} can be simplified to 

k ^ ' u 

where v is any pair partition of {1,2,..., 2n} i— > {1,2,..., n}. As limTv A^~ n (^) = 1/n!, 
the result follows. 



Remarks 

a. Condition (20) is only sufficient to obtain the central limit theorem of Proposition 3. 
The result also follows from weak clustering if the asymptotic state (p^ is permutation 
invariant. Indeed, we could then average first over the time tk p that appears only 
once and use weak clustering to conclude that this average is zero. According to 
Proposition 2, this occurs when (21, G, 4>) is strongly clustering, hence weakly clustering 
and (^oo automatically permutation invariant. 

b. Condition (20) is implied by strong clustering (6. a) because of the equivalence of (6. a) 
and (6.b) and implies weak clustering (5). In fact, with Y := Y — <f>(Y), condition (20) 
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means that 

lim (j)(xY(t)z)= lim (j)(x(-t)YZ(-tj) = 



As far as fluctuations are concerned, permutation invariance (16) implies 

Corollary 1. Let the dynamical system (21,0,0) be weakly clustering. Let the multiple 
average (14) exist and define an asymptotic state (j)^ which is permutation invariant in 
the sense of (16). Then, with . . . ,X^ centred observables of 21, 

(0 r = 2n + 1 

Y m Ux&-x22«) r = 2», 



/ j ^00^2,(1) vv z/(2n) 
ord 

where ord means that the sum is over all ordered pair partitions v = ((oi,/3i), 
(ck2,/#2), • • • ? /9n)) o/ {1, 2, . . . , 2n} , i.e. we choose sites otj < (3j such that v(otj) = 
= j with j running from 1 to n and ct\ < a-i < ■ ■ ■ < a n . 

Proof: An easy consequence of Proposition 3. 



Corollary 2. If the dynamical system (21, 0, 4>) is strongly clustering, then the fluctuations 
Fn(X) of observables X = X* e 21 such that 4>{X) = and <fi(X 2 ) = a 2 tend to Gaussian 
random variables with zero mean and variance a. 

Proof: Since <f>(X) = implies 4>oo(Xj) = 0, for all j G No, we are in fact dealing with 
centred words in Sl^. Because of Proposition 2, we can apply Corollary 1, using the 
notation for pair partitions introduced there, to compute the even moments (the odd ones 
are zero) 

M 2n :=\im ( j> 0O ^F N (X) 2n ^ . 
As in Proposition 2, it follows that 

n 

4>oo (Xv{l)X v (2) ■ ■ ■ X u (2n)^ = Y\_ )) = °" 2n ' 

i=i 

where we have used that 4>oo{x v ^ aj )X v ^.-^j = 0oo^(X 2 )^ = (ft(X 2 ) = a 2 . Since the 

number of pair partitions v : {1, 2, . . . , 2n} i-> {1,2,..., n} is n"= (^T^) = W/ 2 ™; 
we get 

M 2n = (2n-l)lla 2n , 
13 



that is the 2n-th moment of the Gaussian distribution go^(t) with zero mean and variance 
a: g 0i<T (t) = exp(-t 2 /2a 2 ). 



Proposition 3 is sufficient to guarantee that the linear functional $ defined as a point- 
wise limit by (19) is positive on the free *-algebra T generated by all possible fluctuations 
F(X), X G 21, where F(X)* := F(X*). Furthermore, one can look for concrete Hilbert 
space representations of the abstract couples (T , $) . More precisely, one may search for 
a structure determined by creation and annihilation operators a* and a satisfying certain 
commutation relations [14] and for a "ground" state O such that 

(fi,(a + a*) n O) = $(F(X) n ) . 

It is, however, clear that associating definite probability distributions and representations 
to fluctuations very much depends on how well-behaved the even moments in Proposition 3 
are. There are a few cases, Corollary 2 being one of them, where a definite structure 
emerges. We present two of them, but first recall the notion of non-crossing pair partitions. 

Let v be a pair partition of {1,2, ...,2n} into n pairs (otj,(3j), with aj < (3j, such 
that v{czj) = v(/3j) = j, j G {1,2,..., n}. A crossing occurs when 

aj < ctfc < j3j < (3k for some j, k G {1,2, ... ,n} . 

Denoting by c(u) the number of crossings in a pair partition v, v is non-crossing if c(u) = 0. 
If v is a non-crossing pair partition of {1,2,..., 2n}, then its pairs (aj, Pj), j = 1, 2, . . . , n, 
are nested, that is if aj < ak < 13 j for some j, k, then also /3k < f3j. 

• Bosonic Brownian motion 

n 

fe=i 

can be represented in terms of Bosonic creation and annihilation operators, i.e. aa* — 
a*a = 1. 

• Free Brownian motion 

n 

fc=i 

can be represented on the full Fock space in terms of creation and annihilation oper- 
ators satisfying aa* = 1. 

In the case of the Bosonic Brownian motion, the ground state distribution is Gaussian. 
In the case of the Free Brownian motion the ground state distribution is the semicircle law 
and to it only non-crossing pair partitions contribute [7]. Indeed, the semicircle distribution 

7o,iW = 77^V / 4-t 2 , |*| < 2 
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has vanishing odd moments, whereas the even ones are given by the Catalan numbers C n : 
M 2n = { 2 ")/(n + 1) = C n . The C n are defined by the recursion relation 

n 

Co = 1 7 Cn = ^2 Ck-1 Cn-k • (22) 

k=l 

The Catalan numbers count in how many different ways one can partition {1,2,..., 2n} 
in ordered non-crossing pairs {(cti, Pi), (a 2 , P2), • • • , (a n , Pn)}, with ol\ < a 2 < ■ ■ ■ < a n . 

In the following proposition, we exhibit conditions on the dynamical system (21, 0, </>) 
and on the procedure of averaging multi-time correlation functions such that a free statistics 
for time-asymptotic fluctuations arises. We start by observing that the explicit rule (14) 
for computing averages has a natural product structure: if t 1— > f\ (t^ , t i2 , . . . , t im ) and 
t 1— > f2(tj 1 ,tj 2 , . . . ,tj n ) are uniform limits of multi-time correlation functions, then 

{z 1 ,z 2 ,...,z m }n{j 1 ,j 2 ,..., Jn } = =}► Avg(/ 1 / 2 )=Avg(/ 1 ) Avg(/ 2 ) . (23) 

Notice, however, that an average is generally not specified if we impose the product prop- 
erty and specify all single-time averages Avg. Indeed, e.g. a bounded function / over Z 2 is 
not necessarily the uniform limit of linear combinations of products of bounded functions 
ti,t 2 h-> h(ti)g(t 2 ) over Z. 

Proposition 4. Let (21, G, 4>) satisfy condition (20) and let the multiple average (14) 
exist and define an asymptotic state (j)^ on SIoq. Let us also assume that for any time 
independent choice of observables A and C and centred observables X , Y and B 

Avg(t ^ (j>(AX(t)BY(t)Cy^ = . (24) 
Then, if XW,...,X( 2n ) are centred observables in 21, 

Ti 

^j^{f n {x^)f n {x^)...f n {x^)) = ^£ <W) UK x{ak)xiPk) ) ' 

v k=l 

where the sum extends over all pair partitions v = ((cti, Pi), (a 2 , P2), ■ ■ ■ , (a n , Pn)) of 
{1,2,..., 2n}. If the asymptotic state <poo is also permutation invariant, then 

^^(F N (X^)F N (X^)...F^ , 

v, ord k=l 

where the sum now extends over all ordered pair partitions v of {1, 2, ... , 2n}. 

Proof: As the asymptotic state </>oo is defined according to (14), we know it to be strongly 
compatible and endowed with the product structure (23). It need not be, however, per- 
mutation invariant. In any case, on the basis of Proposition 3, we only have to consider 
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even moments of order 2n. More precisely, we can restrict to expectations of the form 
(poo (^(i) • • • -X^(2n)) ' wnere ^ is a pair partition of {1,2,..., 2n}. In order to compute 
such an expectation, we must perform successive time-averages as in (14). We begin with 
averaging with respect to t\ 

Either the time t\ appears in two consecutive words X^ and X^ +1 \ in which case we 
can use weak clustering to obtain 

Avg(f! ~ ^{x^\t v{1) )X^\t v{2) ) ■ ■ ■ X^ n \t v{2n) ))) 

= <p(x^x^) • • ■x^\t v{j _ 1) )x^\t v{3+2) ) ■ ~x<M(t v{2n) j)) . 

Or, the two words at time t\ are separated by one or more other words at different times: 
cj)(w X {ai \ti) w'X^ih) w'^j. In this case we use condition (24) as follows. Notice that 
w' is not necessarily centred and so we write 

w ' = (w' - (p(w')l) + (p(w')l . 

Then, remembering that the time t\ does not appear in w, w', w", condition (24) and 
weak clustering imply 

Avg(ti i-> <j>{w X^\t x ) w'X^Xh) «;")) 

= <j)(w') Awg(t 1 ^(l>(w(x^X^yt 1 )w ,,> j > j 

= (j){X {a ^X { ^) <j)(w') <j)(ww") . (25) 

Consider now the average with respect to t 2 . This is either a similar average as that with 
respect to t\ or else the time t 2 appears separately in each of the factors <j>(w') and 4>(ww") 
in (25). If so, the average with respect to t 2 vanishes because of weak clustering and the 
centredness of the X's. Repeating this argument for all the time averages, we see that 
only nested pairs i.e. non-crossing pair partitions, contribute whence the result. 



Proposition 4 and permutation invariance of the asymptotic state provide us with 
sufficient conditions on the correlation functions to yield fluctuations which are free random 
variables and hence semicircular ly distributed. 

Corollary 3. Let the dynamical system (21,0,0) satisfy conditions (20) and (24) and 
the multiple average (14) exist and define a permutation invariant asymptotic state (j)^ as 
in (16). Then, the fluctuations Fn(X) of observables X = X* e 21 such that 4>{X) = 
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and 4>(X 2 ) = a 2 tend to semicircularly distributed random variables with zero mean and 
variance a. 

Proof: We can apply the previous proposition to compute the even moments (the odd 
ones are zero) 

M 2n := lim^FjvpT) 2 ™) = C n a 2n , 

where C n are the Catalan numbers (22). On the other hand, the latter moments are the 
non-zero moments of the semicircle distribution 70, CT (t) = ^/4a 2 — t 2 /2jva 2 on \t\ < 2a. 

m 

Remarks 

a. Referring to the types of clustering we presented in Section 2, clustering in the 
mean (4) would not be sufficient to yield the result of Proposition 4, while strong 
clustering (6. a) is incompatible with condition (24) and with freeness, as already ob- 
served in Section 3. In the proof of Proposition 4, weak clustering (5) has been used 
in an essential way. Indeed, suppose that we assume instead of weak clustering only 
clustering in the mean. In the process of centring the observable «/, sandwiched be- 
tween two equal time observables X^ k \tk) and X^ k \tk), some other time tt might 
appear once in both w' and in to". Therefore, we should consider an average of a 
product 1 1— > f(t)g(t) which is generally unequal to the product of the averages 

Avg(t » f(t)g(tj) ^ Avg(t » /(*)) Avg(t » . 

b. Conditions (20) and (24) are independent. Obviously, the latter does not imply the 
first one. As a counterexample to (20) =>- (24), consider the algebra 21 generated by 
linear combinations of the unitary operators W(p), p = (^1,^2) £ Z 2 , such that 

W(p)* = W(-p) , W(p)W(q) = e l7T0a{p ^W(p + q) , 

where 9 G [0,1) and a(p,q) = p\q 2 — P2Qi- By equipping 21 with the tracial state 
4>(W(p)) = 5 P; o and the automorphisms 0*(W(p) = W(T*p), where T is the typical 
2x2 matrix of the hyperbolic dynamics on the two-dimensional torus, one obtains 
the class of quantized cat-maps studied in [15]. A thorough examination of such 
models in the framework presently developed will be the subject of a forthcoming 
paper [16]. Here, it is sufficient to mention that there are values of 9 for which 
lim t 9a{T t p, q) mod 1 = /3A(p, q) with j3 and A(p, q) two not necessarily zero quan- 
tities. 

One calculates (p(w(T^p)W(T t2 ci)W(-p)W(-T t2 q))= e 2 ™ e ^ Tt2q ' p )5 T t 1P) p. Thus, 
Avg(t ^ <p(w(p)W(T t q)W(-p)W(-T t q)))=e 27Tlf3A ^' p \ whereas 
lim tl it2 ^ oo 0(w(T*ip)W(r t »q)W(-p)W(-T*»q))=O. 
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5 Examples 

We implement in this section the constructions of above for a few simple model systems. 
More complex dynamical systems, with richer behaviour of correlation functions, will be 
considered in a forthcoming paper [16]. As Proposition 2 and Corollary 1 deal with the 
asymptotic state and the fluctuations of strongly clustering dynamical systems, we shall 
only consider here more extremely non-commuting dynamical systems. 

As a first model, we consider a probability space (X, fx) equipped with an automor- 
phism 9 i.e. a measurable transformation 9 of X, with measurable inverse 9~ l , such that 
fi = fio9 = fio 6> _1 . We will assume that 9 is mixing 

hm fu(fgo6 n )=iu(f)iu(g) for all /, g e C 1 (X, . 

n— >oo 

In the Koopman description of the classical dynamical system (X, fj,, 6>), one considers the 
Hilbert space ft := C 2 (X, fx) and the unitary single step evolution Uf := f o 9 of the "wave 
functions" /. The expectation fj,(f) of an C°° function on X is recovered as 

/*(/) = (1,M/1> . 

In this formula 1 denotes the constant function 1 on X and Mf the multiplication oper- 
ator on f) by the function /. We now consider a "non-canonical" quantization whereby 
the quantum evolution of "wave-functions" is exactly the classical one [17]. In such a 
description, one may choose for 21 the algebra of compact operators CI + &(Sj) to which 
a unit has been added. The dynamics is determined by the usual Heisenberg picture: 
Q(X) := UXU* and the reference state <fi is the vector state defined by 1. Using the 
assumed mixing property of (X, li, 9), one readily verifies that whenever X e &($)) and w 
and w' are products in < 2i OCl of copies of compact operators, such that both the rightmost 
letter of w and the leftmost of w' are different from j, then 

^(wXjw') = <j>{X) (PUw) tooiw') , (26) 
even when w or w' contain letters pertaining to the index j. As a consequence, 

n 
3 = 1 

whenever X^ e &(S)) and (p^ is permutation invariant. For a compact X e &(S)) we will 
denote by X the element X — <j>(X)l obtained by centring X. In spite of good ergodic 
properties of (X, 9, fi), its "quantization" (21, 6, (j>) is not weakly clustering as can be seen 
from 

hm cj>(XQ t (Y)X) = <P(X) 2 0(Y) ? <p(X 2 ) cj>{Y) , 
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X and Y compact. Considering fluctuations is not meaningful in such a case. In fact, 
expectations of third order moments of fluctuations F N {X) diverge with N as can be seen 
from 



'N VN 
with X :=X -(p(X)l. 

As a second example, we consider a Fermionic dynamical system whose observables 
are given by a CAR-algebra %l(Sj) over a single-particle space ft. 21(f)) is generated by an 
identity 1 and by creation and annihilation fields {a*(<p) \ <p e $)} and {a (if) \ <p G $)} 
subject to the relations: 

y> i— > a*(y>) is C-linear and 

{a((p),a(ip)} = and {a(p), a*(ip)} = ((p, ip) . 

The parity automorphism n on 21(f)) is determined by 7r(a*(<^)) := — a* ((f) and a *- 
automorphism is said to be even if o n = n o 0. We assume furthermore that the 
dynamics determined by such an even is asymptotically "Abelian" in the sense that 



lim ||{a*(</?),0((a # 







where denotes either a or a*. It is well-known that a reference state 0, invariant under 
an even, asymptotically Abelian 0, is automatically even i.e. o tt = 0. Finally, we shall 
assume that the dynamical system is multi-clustering 



k 

lim <f>(xW(t 1 )xW(t 2 )---XW(t k j) = TT<K* (j) ) 

tj—tj + l—HX V / - LJ - 

j = l 



with ti > t 2 > • • ■ > tk- The asymptotics of multi-correlation functions can be computed 
along the same lines as that for the strongly clustering case and we obtain when each of 
the X^ is either even or odd 

...*<& ) - n^oL * w ) • (^) 

P =i p 

I p is, as in Proposition 2, the set of all indices j such that = p and e is either 

1 or —1 according to whether an even or odd permutation is needed to permute the 
odd XW appearing in X^^X^^ • ■ • X^^ into the order in which they appear in the 

product at the right-hand side of (27). This amounts to saying that (p^ is the Chevalley 
product of a countable number of copies of 4> on the twisted tensor product 21(© N(I 5]) of 
No copies of 21(f)) with itself. The fluctuations of creation and annihilation fields are now 
straightforwardly computed, yielding 

lim^(/) Q0 ^F N (a*(^ 1 ))F N (a*(^2)) ■ ■ ■ F N (a*(^ n ))^j = QF (a # (v?i)a # (</? 2 ) • • • a # (</?n)) • 



(f>QF is the "quasi-free projection" of the state 0, namely the quasi-free state on 21(53) 
determined by the covariance 

(<£,V0 i-> </>(a # (^)a # (^)) . 
We recover hereby the Fermionic central limit theorem of [18]. 

In its most basic form, the "free shift" is a quantum shift 

6(e;) := e i+ i, z G Z 

on a set of generators {e^ | z G Z} that satisfy the relations 

e* = and ef = 1, i G Z . 

The algebra 21 of observables is the universal C*-algebra generated by 1 and the e,. The 
finite linear combinations of monomials of the type ej 2 • • • ei m with ^ i^+i form a dense 
*-subalgebra 2lo of 21. The product of two monomials determined by ordered index sets 
{z'i, 12, . . . , z' m } and {ji, j2, • • • , in} is the monomial corresponding to the index set obtained 
by first concatenating {z'i, i 2 , . . . , i m } and {ji,j2, • • • , jn} and then omitting those indices 
that appear twice in subsequent positions. The identity corresponds to the monomial with 
empty index set and the adjoint of a monomial is the monomial with reversed index set. 
If there are no preferred observables to single out apart from the identity, a meaningful 
reference state is the state on 21 

(p(e tl e t2 ■■■e im ) = if m > and 0(1) = 1 . (28) 

The dynamical system (21, 0, (p) is weakly, but not strongly, clustering. It is quite straight- 
forward to compute the asymptotic state (j)^: it is the free product *j 6 N </> of copies of (f>. 
As each centred element in 2lo is a finite linear combination of monomials with non-trivial 
dependence set, it suffices to show that for n = 1,2, ... the expectation of an element 
^(1)^(2) ' ' ' ™{n) m tne state 0°° vanishes, where each X^ is a non-trivial monomial 
and consecutive v(k) are different. Clearly, when all differences |t& — tg\ for k ^ I appear- 
ing in the list of ^(j)'s become sufficiently large, there is no possible simplification in the 
monomial X^(t u ^)X^(t u ^ 2 )) • ■ ■ X^ (t„( n )) due to the rule ej = 1 and because of (28) 

cj>{X^{t <x) )XV){t v{2) ) ■ ■■X^){t v{n) )) = and so 0oo(*$)*i(2) ' ' '4(n)) = to °- The 
state is permutation invariant and satisfies even a strengthened Condition (24), where 
the average is replaced by a limit. Temporal fluctuations of centred self-adjoint observables 
are therefore, by Proposition 4, semicircularly distributed. 
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